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1.  INTRODUCTION 

Stellar  scintillation,  defined  here  generally  as  the  inten¬ 
sity  fluctuations  of  a  plane  wave  after  propagating  from 
space  to  Earth,  is  typically  quantified  in  terms  of  the  ir- 
radiance  variance  at  a  point  on  a  telescope  aperture, 

-  {r^)l{I)^  -  1.^  In  the  regime  of  weak  turbulence  (i.e., 
(.7/^  <  1),  scintillation  is  well  understood  and  described  by 
Rytov’s  theory,  where  c7/^  just  depends  on  the  total 
altitude-weighted  turbulence  along  the  propagation 
path^: 


ar  =  (Tr 


exp 


2,25^2 


-  1.  (1) 


Here  the  refractive  index  structure  constant  C„^(/i)  is  a 
measure  of  turbulence  strength  at  altitude  h,  <f)  is  the  ze¬ 
nith  angle,  and  k  is  the  optical  wave  vector’s  magnitude. 
However,  in  the  regime  of  strong  turbulence  it  is  well 
known  that  the  relationship  betw'een  rr/“  and  breaks 
down.  Specifically,  for  (r/^  >  1  the  irradiance  variance 
saturates,  that  is,  it  reaches  a  peak  (refen-ed  to  as  the 
strong-focusing  regime)  and  then  slowly  decreases  with 
increasing  total  turbulence.^ 

Over  the  past  decade  considerable  progress  has  been 
made  in  understanding  the  saturation  of  scintillation,  in 
particular  for  horizontal  propagation  through  homoge¬ 
neous  turbulence.  For  example,  the  role  of  the  inner 
scale  has  been  elucidated'^  and  has  been  found  to  have  a 
greater  significance  for  spherical  waves  than  for  plane 
waves.'’’  Additionally,  the  probability  density  function  of 
irradiance  has  been  numerically  evaluated,  and  its  devia¬ 
tions  from  a  log-normal  distribution  have  been  studied.^ 
The  situation  is  quite  different  with  regard  to  stellar  scin¬ 
tillation,  however,  where  very  little  is  known  about  satu¬ 
ration.  In  large  part  this  is  because  the  irradiance  vari¬ 


ance  typically  associated  with  stellar  scintillation  is  small 
enough  that  Rytov’s  theory  applies.'  Nevertheless,  de¬ 
lations  from  Rytov’s  theory  can  occur,  especially  at  large 
zenith  angles,^  and  so  it  is  reasonable  to  question  how  the 
saturation  of  stellar  scintillation  might  differ  from  that 
associated  with  propagation  through  homogeneous  turbu¬ 
lence,  In  particular,  since  expression  (1)  is  no  longer 
valid  in  the  case  of  saturation,  one  might  wonder  whether 
a/-  still  depends  on  just  the  integrated  turbulence  along 
the  entire  propagation  path  or  whether  certain  atmo¬ 
spheric  regions  play  different  roles  in  producing  scintilla¬ 
tion. 

In  the  following  sections,  we  consider  this  question  by 
numerically  simulating  optical  propagation  through  two 
distinct  profiles  of  atmospheric  turbulence.  For  one  pro¬ 
file  scintillation  is  dominated  by  turbulence  in  the  tropo- 
pause,  while  for  the  other  both  tropopause  and  low- 
altitude  tropospheric  turbulence  contribute  to 
scintillation.  In  Section  2  we  briefly  discuss  the  numeri¬ 
cal  simulation  procedure  and  our  vertical  profiles  of  atmo¬ 
spheric  turbulence.  The  results  are  presented  in  Section 
3,  where  the  combination  of  high-  and  low-altitude  turbu¬ 
lence  is  shown  to  be  more  effective  in  producing  scintilla¬ 
tion  in  the  strong-focusing  regime  than  high-altitude  tur¬ 
bulence  alone.  In  Section  4  w^e  explain  this  result  by 
considering  the  relationship  between  a  propagated  field’s 
lateral  coherence  length  and  its  Fresnel  length  for  the  two 
different  atmospheric  profiles. 


2.  VERTICAL  PROFILES  OF  ATMOSPHERIC 
TURBULENCE  AND  THE  NUMERICAL 
PROCEDURE 

To  describe  the  atmosphere’s  vertical  turbulence  profile, 
we  employ  a  slightly  modified  version  of  the  Hufnagel- 
Valley  (HV)  modeP: 
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C,,Hh)  =  r{8.15  X  10”^^i;//i^®exp[-(/i/1000)]  +  2.7 

X  10“^^  exp[ -(/i/lSOO)]  +  A  exp[-(A/100)]}. 

(2) 

Here  is  a  parameter  in  meters  per  second  that  controls 
high-altitnde  turbulence,  and  A  accounts  for  boundary- 
layer  atmospheric  turbulence.  (In  the  literature,  is 
sometimes  defined  as  the  rms  wind  speed  in  the  5-20-km 
altitude  range.^)  F  is  a  factor  that  scales  a  particular 
profile  to  higher  and  lower  levels  of  turbulence  and  simul¬ 
taneously  could  account  for  the  sec^^^(  6)  zenith-angle 
factor  in  expression  (1). 

A  particular  profile  of  turbulence  that  is  often  used  in 
studies  of  optical  propagation  through  the  atmosphere  is 
the  HV-21  profile;  =  21  m/s  and  A  =  1.7 

X  10” Here,  we  examine  this  profile's  influence 
on  the  saturation  of  scintillation  by  increasing  F  for  a 
1-^m  reference  wavelength.  In  this  way,  we  keep  the 
relative  shape  of  the  profile  constant,  while  increasing  the 
level  of  integrated  turbulence  (i.e.,  The  altitude- 

weighted  contribution  of  the  HV-21  profile  to  stellar  scin¬ 
tillation  [i.e.,  C^\h)h^^]  is  shown  in  Fig.  1  for  the  specific 
case  of  =  0.5,  and  it  is  apparent  that  both  high- 
altitude  (-11  km)  and  low-altitude  (-1.4  km)  turbulence 
contribute  to  scintillation.^  For  comparative  purposes, 
we  have  also  considered  a  profile  that  we  refer  to  as  the 
high-altitude  (HA)  profile.  For  the  HA  profile,  F  =  1, 
and  is  varied  to  produce  different  levels  of  turbulence; 
A  is  somewhat  arbitrarily  set  to  1.0  X  10"^^ 
Though  the  HA  profile  may  appear  somewhat  artificial,  as 
will  be  seen  it  provides  a  good  countertype  for  the  more 
realistic  HV-21  profile.  Basically,  the  idea  behind  the 
HA  profile  is  to  have  a  point  of  comparison  that  will  allow 
us  to  better  understand  those  features  of  a  turbulence 
profile  that  play  a  significant  role  in  stellar  scintillation. 
The  HA’s  altitude-weighted  contribution  to  stellar  scintil¬ 
lation  is  also  shown  in  Fig,  1  for  the  case  of  =  0,5, 
and  it  is  clear  that  scintillation  in  this  case  is  influenced 
primarily  by  turbulence  in  the  tropopause. 

Propagation  through  the  turbulent  atmosphere  is  stud¬ 
ied  numerically  according  to  standard  methodology. 
Briefly,  we  simulate  three-dimensional  wave  propagation 


Fig.  1.  versus  altitude  for  the  Hufnagel- Valley 

model  of  atmospheric  turbulence.  For  the  HV-21  profile,  A 
=  1.7  X  lO"  m  i\  =  21  m/s,  and  the  strength  of  turbu¬ 
lence  is  varied  by  changing  F  (for  the  graph  F  =  3.85).  For  the 
high-altitudc  profile,  labeled  HA,  A  =  1.0  x  10  ^^m  F  1, 
and  the  strength  of  turbulence  is  varied  by  changing  the  rms 
wind  speed  (for  the  graph  u,  =  47  m/s).  Each  profile  shown  in 
the  figure  yields  (7^'^  ~  0.5. 


by  placing  a  number,  ,  of  thin  random-phase  screens  at 
specific  altitudes,  h,.  Prior  to  impinging  on  the  ith 
phase  screen,  the  optical  field  E(p,hp  is  given  by 

h~)  =  EJp,  h~)exp[l(b(p,  /if)],  (3) 

where  p  is  a  position  vector  in  a  plane  perpendicular  to 
the  propagation  direction  (for  the  computations,  the  ze¬ 
nith  angle  is  zero),  and  where  and  0  are  real  scalar 
fields.  The  phase  screen  imparts  an  additional  random 
phase  to  this  wave  d  (where  the  statistical  properties  of  0 
are  determined  by  the  characteristics  of  turbulence  at  hi), 
so  that  after  propagation  through  the  phase  screen,  the 
field  becomes  E{p,  /if)  =  E^{p,  /if )exp{i[d<p,  /if) 
-  //(p, /if)]}.  Taking  the  two-dimensional  Fourier 
transform  of  E(p,  /if),  the  field  is  propagated  to  the  next 
screen  under  the  Fresnel  approximation,^^  and  after  an 
inverse  Fourier  transform  this  becomes  E{p,  /if^i).  We 
start  from  an  initial  monochromatic  plane  wave  at  an  al¬ 
titude  of  20  km,  and  the  wave  is  propagated  to  a  telescope 
aperture  on  the  Earth's  surface.  More  details  on  the  nu¬ 
merical  simulation  may  be  found  in  Appendix  A, 

A  single  realization  consisted  of  plane-wave  propaga¬ 
tion  (\  =  1  pm)  from  the  lower  stratosphere  to  the  Earth 
through  four  512  X  512  pixel  phase  screens  of  equal 
length.  Depending  on  the  simulation  parameters,  the 
phase-screen  lengths  ranged  from  170  to  210  cm.^^  For 
each  simulation,  we  recorded  the  field  intensity  at  the 
central  pixel  of  a  25-cm-diameter  telescope  aperture, 
and  we  determined  the  product  of  the  central  pixel’s  log 
amplitude  with  the  log  amplitude  of  all  other  aperture 
pixels.  Additionally,  we  determined  the  aperture- 
averaged  intensity  of  the  propagated  field  for  each  simu¬ 
lation.  A  single  computational  run  involved  3000  real¬ 
izations  of  atmospheric  propagation,  and  four  to  seven 
runs  were  used  to  generate  pooled  estimates^^  of  irradi- 
ance  variance.  The  log-amplitude  correlation  functions 
were  evaluated  from  a  single  computational  run  of  15,000 
realizations.  Though  one  could  attempt  to  calculate 
using  fewer  simulations  (for  example,  by  treating  each 
pixel  from  a  single  realization  as  a  separate  estimate  of 
the  intensity  fluctuations),  we  chose  a  different  approach, 
since  the  intensity  fluctuations  of  neighboring  pixels  are 
correlated.  Calculations  were  performed  on  a  dual¬ 
processor  600-MHz  DEC  Alpha  Server  4000,  and  a  single 
realization  required  —18  s. 

As  evidence  of  the  numerical  procedure's  validity,  Fig, 
2  shows  the  log-amplitude  correlation  function  [i.e., 
B^(p)/B^{0)]  for  the  case  of  HV-21  turbulence  with  F 
=  1.  Since  the  turbulence  is  relatively  weak  in  this  case 
(i.e.,  =  0.11),  the  numerical  results  may  be  compared 

with  Rytov  theory,  where 


ri!  r 

BJp)  =  C,;\h)  (0.033) 

Jo  Jo 


o(xp)K* 


X 


K^h  \ 


KdK 


dh. 


(4) 


In  writing  Eq.  (4),  we  have  restricted  ourselves  to  Kol¬ 
mogorov  turbulence.  The  solid  curve  in  Fig.  2  coire- 
sponds  to  the  prediction  of  Eq.  (4),  and  as  one  would  ex- 
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Fig.  2.  Circles,  log- amplitude  correlation  function  for  the  HV-21 
model  with  T  =  1  (i.e.,  weak  turbulence):  -  11.8  cm,  and 

=  0.11.  Solid  curve,  prediction  from  Rytov’s  theory  with 
Eq.  (4). 

pect,  the  agreement  between  Rytov  theory  and  the 
numerical  results  is  quite  good. 

3.  RESULTS 

Our  principal  results  are  shown  in  Figs.  3  and  4,  where 
(Tj^  and  the  aperture-averaged  irradiance  variance 
ai*e  shown,  respectively,  as  functions  of  the  Rytov  stan¬ 
dard  deviation  an .  Open  circles  in  the  figures  corre¬ 
spond  to  the  HV-21  profile,  solid  circles  correspond  to  the 
HA  profile,  and  error  bars  indicate  the  variance’s  stan¬ 
dard  error  as  estimated  from  the  four  to  seven  computa¬ 
tional  runs.  (As  a  point  of  reference,  we  note  that  a^ 
=  2  for  the  HV-21  profile  can  be  attained  by  attributing 
r  solely  to  a  77°  zenith-angle  factor.)  As  a  test,  we  per¬ 
formed  several  runs  for  a  numerical  model  with  five  phase 
screens  and  found  no  statistically  significant  difference  in 
the  results.  The  dashed  curv^e  in  Fig.  3  con-esponds  to 
solid  curves  are  fits  of  the  data  to  the 

empirical  formula: 

<Ti'  =  -  exp[-(aK^/(r„-}]].  (5) 

Here  is  a  parameter  that  indicates  both  the  level  of 
saturation  in  the  strong-focusing  regime  and  the  Rytov 
variance  where  saturation  sets  in.  It  should  be  noted 
that  Eq.  (5)  is  not  based  on  any  underlying  theory  of  scin¬ 
tillation;  rather,  it  is  simply  a  convenient  analytical  form 
for  discussion  of  the  numerical  results. 

Regarding  Fig.  3,  it  is  apparent  that  the  two  profiles 
display  different  saturation  characteristics.  Specifically, 
for  the  HA  profile  in  the  strong-focusing  regime  is 
smaller  than  for  the  HV-21  profile.  In  terms  of  Eq. 
(5)  this  may  be  summarized  as  (7^,^(HA)  =  1.6  and 
c7^‘ (HV-21)  =  2.5.  Regarding  Fig.  4,  we  have  the  alter¬ 
native  finding  that  (a/^)  in  the  saturation  regime  is 
smaller  for  the  HV-21  profile  and  larger  for  the  HA  pro¬ 
file.  These  apparently  contradictory  results  may  be  rec¬ 
onciled  if  the  HV-21  profile  produces  a  log-amplitude  cor¬ 
relation  function  with  a  relatively  short  correlation 
length.  A  short  correlation  length  implies  not  only  sig¬ 
nificant  amplitude  fluctuations  at  a  point  from  instant  to 
instant  (i.e.,  large  a/^  for  the  point  intensity  from  realiza¬ 
tion  to  realization),  but  also  a  greater  ability  for  the  tele¬ 
scope  aperture  to  average  out  amplitude  fluctuations  at 
any  given  instant  (i.e.,  small  (af)).  Figure  5  shows  the 
log-amplitude  correlation  functions  for  the  two  profiles  in 


the  case  of  a^  =  4.8,  and  it  is  clear  that  the  HV-21  profile 
has  the  shorter  correlation  length.  Consequently,  Figs. 
3-5  are  self-consistent  and  result  in  two  immediate  con¬ 
clusions.  First,  the  atmosphere’s  vertical  turbulence  pro¬ 
file  has  an  important  influence  on  the  saturation  of  stellar 
scintillation.  Second,  the  combination  of  high-  and  low- 
altitude  turbulence  is  more  effective  in  producing  scintil¬ 
lation  in  the  strong-focusing  regime  than  high-altitude 


Fig,  3.  Phase-screen  calculation  of  irradiance  variance  ver¬ 
sus  the  Rytov  standard  deviation  (Xn  .  Open  circles,  HV-21 
model;  solid  circles,  HA  model.  Dashed  curv^e,  Rytov  theory  fi.e., 
expression  (1)];  solid  cur\^es,  Eq.  (5)  with  =  1.6  and 
=  2.5  for  the  HA  and  the  HV-21  profiles,  re.spectively. 


Fig.  4.  Phase-screen  calculation  of  the  aperture-averaged  irra¬ 
diance  variance  (cr/^)  versus  the  Rytov  standard  deviation  , 
Open  circles,  HV-21  model;  solid  circles,  HA  model.  The  tele¬ 
scope  diameter  in  the  simulations  was  25  cm. 


Fig.  5.  Log-amplitude  correlation  function  B p)/B ^.{0)  in  the 
strong-focusing  regime  (i.e.,  cr^  -  4.8):  open  circles,  HV-21  pro¬ 
file;  solid  circles,  HA  profile. 
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turbulence  alone.  We  present  a  physical  explanation  for 
these  conclusions  in  the  next  section. 


4.  DISCUSSION 

As  a  way  to  explain  scintillation  qualitatively/^  turbulent 
eddies  are  often  imagined  as  extremely  long  focal  length 
lenses  that  randomly  focus  and  defocus  a  coherent  wave. 
For  eddies  of  size  I  at  an  altitude  /i,  the  focal  length  varies 
like  Moreover,  in  the  regime  of  weak  turbu¬ 

lence,  where  the  effects  of  the  individual  eddies  are  inde¬ 
pendent,  the  variance  of  the  fluctuations  scales  like  the 
number, of  eddies  along  the  propagation  path.  Thus,  not 
only  do  small  eddies  produce  the  greatest  intensity  fluc¬ 
tuations  at  the  point  of  observation,  but  also  there  are 
more  of  them  along  the  propagation  path.  Taken  to¬ 
gether,  these  statements  imply  that  the  smallest-scale  ed¬ 
dies  dominate  the  production  of  scintillation.  Of  course, 
the  smallest  eddy  that  can  efficiently  focus  or  defocus  the 
light  in  this  fashion  is  constrained  by  diffraction,  which 
limits  the  smallest-scale  eddies  of  relevance  for  scintilla¬ 
tion  to  the  Fresnel  length. 

While  the  above  description  of  scintillation  is  appropri¬ 
ate  in  the  regime  of  weak  turbulence,  and  in  fact  provides 
a  semi  quantitative  rationale  for  the  Rytov  variance’s 
form,  it  lacks  an  important  subtlety.^*^’^®  In  the  regime  of 
strong  turbulence  a  wave  is  perturbed  by  its  propagation 
through  the  turbulent  atmosphere  before  interacting  with 
a  Fresnel-scale  eddy,  with  the  consequence  that  this  eddy 
must  now  focus  or  defocus  a  partially  coherent  field.  Be¬ 
cause  the  lateral  coherence  length  p„  of  a  propagated  field 
decreases  as  the  path-integrated  turbulence  increases,  a 
point  is  eventually  reached  at  which  the  coherence  length 
becomes  smaller  than  the  Fresnel  length.  At  this  level  of 
turbulence,  the  Fresnel  scale  eddies  lose  their  ability  to 
effectively  focus  or  defocus  the  wave.  Consequently,  fur¬ 
ther  increases  in  turbulence  do  not  produce  much  in¬ 
crease  in  scintillation,  and  the  irradiance  variance  satu¬ 
rates. 

In  the  present  problem,  as  suggested  by  Eq.  (1)  and  Fig. 
1,  the  most  significant  Fresnel  length  for  scintillation  is 
associated  with  the  tropopause  for  both  the  HA  and  the 
HV-21  profiles.  Consequently,  scintillation  should  satu¬ 
rate  when  the  field’s  coherence  length  after  propagation 
from  space  to  the  tropopause  becomes  smaller  than  the 
tropopause  Fresnel  length  [i.e.,  Pf,(/i)  <  =  10  cm,  as¬ 

suming  a  l-/xm  wavelength  and  a  10-km  tropopause  alti¬ 
tude].  Clearly,  the  HA  profile  has  more  stratospheric 
turbulence  than  the  HV-21  profile.  Thus,  for  the  same 
degree  of  weighted,  path-integrated  turbulence  from 
space  to  Earth  (i.e.,  rr/j),  the  HA  profile  will  always  pro¬ 
duce  a  smaller  coherence  length  at  the  tropopause,  and 
this  in  turn  implies  that  the  HA  profile  will  saturate  at 
lower  values  of  (tj^  . 

This  argument  can  be  placed  on  a  more  quantitative 
footing  by  writing  the  coherence  length  at  an  altitude  h 
for  a  wave  propagating  from  space  to  Earth  as 


1.45/e“|  C,^{/l')d/^' 

h 


--3/0 


(6) 


Scintillation  parameter,  f 


Fig.  6.  Rytov  variance  parameterized  by  the  scintillation  pa¬ 
rameter  f  as  discussed  in  the  text.  Open  circles,  HV-21  profile 
with  /“averaged  over  the  tropopause;  solid  circles,  HA  profile  with 
f  averaged  over  the  tropopause;  dashed  curve,  HV-21  profile  with 
f  averaged  over  the  lower  troposphere.  For  the  tropopause- 
averaged  curves,  the  HV-21  and  the  HA  profile  Rytov  variances 
at  /*  =  1  are  in  reasonable  agreement  with  their  respective  cr„“ 
values. 


Forming  the  ratio  of  p„(/i)  with  the  Fresnel  length  and 
then  averaging  over  the  tropopause  region  defines  a 
scintillation-saturation  parameter  f\ 


/•- 


r/>2 

Po(^) 

J/., 

\j\h 

(7) 


where  h  ^  and  h  i  correspond  to  the  top  and  bottom  of  the 
tropopause,  respectively.  When  f  is  less  than  unity  for 
some  vertical  turbulence  profile,  the  tropopause  coher¬ 
ence  length  is  less  than  the  Fresnel  length,  and  scintilla¬ 
tion  should  saturate.  Parametrizing  the  Rytov  variance 
by  f  and  writing  this  as  cr^if),  we  then  expect  o’/  to  satu¬ 
rate  at  Cr;j“(l). 

Figure  6  shows  versus  f  for  the  HV-21  profile 

(open  circles)  and  the  HA  profile  (solid  circles).  As  is 
clear  from  the  figure,  (r^^(l)  for  the  HA  profile  is  smaller 
than  for  the  HV-21  profile,  consistent  with  the  ob¬ 

servation  that  cT/^  saturates  at  lower  values  for  the  HA 
profile.  Moreover,  we  obtain  crj^^d)  ==  2.5  for  the  HV-21 
profile  and  =  1.3  for  the  HA  profile,  which  are 

values  that  are  consistent  with  HV-21)  and  rro^(HA), 
respectively.  For  completeness,  the  dashed  curve  in  the 
figure  shows  versus  /*for  the  HV-21  profile,  but  now 
with  f  averaged  over  the  lower  troposphere  (i.e.,  1-2.5 
km).  Comparison  of  the  two  HV-21  curves  illustrates, 
that  saturation  is  primarily  associated  with  the  tropo¬ 
pause  for  this  profile,  even  though  Fig.  1  suggests  that  the 
lower  troposphere  plays  a  significant  scintillation  role  in 
the  regime  of  weak  turbulence. 


5.  SUMMARY 

In  this  work  we  investigated  the  relationship  between 
stellar  scintillation  in  the  strong-focusing  regime  and  the 
atmosphere’s  vertical  turbulence  profile.  Specifically,  us¬ 
ing  the  Hufhagel-Valley  (HV-21)  model  of  atmospheric 
turbulence,  we  compared  the  irradiance  variance  for  the 
standard  HV-21  turbulence  profile  with  values  ob¬ 
tained  from  what  we  termed  a  high-altitude  (i.e.,  HA)  pro¬ 
file.  In  the  HV-21  profile,  scintillation  arises  from  both 
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low-  and  high-altitude  turbulence,  while  in  the  HA  profile 
scintillation  is  dominated  by  high-altitude  turbulence 
alone.  Comparing  the  values  for  the  two  profiles, 
which  were  obtained  by  numerical  simulation,  we  found 
that  in  the  strong-focusing  regime  the  vertical  profile  of 
turbulence  plays  an  important  role  in  scintillation. 
Moreover,  we  found  that  the  combination  of  high-  and 
low-altitude  turbulence  was  more  effective  in  producing 
scintillation  in  the  strong-focusing  regime  than  high- 
altitude  turbulence  alone.  These  results  were  explained 
as  a  consequence  of  an  optical  wave’s  decoherence,  rela¬ 
tive  to  the  tropopause  Fresnel  length,  as  it  propagates 
through  the  lower  stratosphere. 


APPENDIX  A 

In  the  numerical  simulation  of  optical  propagation,  the 
placement  and  strength  of  the  phase  screens  is  deter¬ 
mined  by  a  Gaussian  quadrature  procedure.^^  Basically, 
since  we  need  to  replace  the  continuous  vertical  turbu¬ 
lence  profile  with  a  set  of  discrete  turbulent  layers,  we  ap¬ 
proximate  the  integral  of  with  a  discrete  sum  of 
phase  screens  by  using  Gauss-Legendre  quadrature 

fH  ro.skm  CH 

=  C„Hh)dh  +  C„hh)dJi 

Jo  Jo  J  0.3  kin 

=  a,C„^{hO  +  2  aiC.Hhi).  (Al) 
1=2 

Here,  H  corresponds  to  a  somewhat  arbitrary  20  km 
height  of  the  atmosphere,  hi  equals  0.15  km,  and  the  hi 
and  ai  in  the  second  term  on  the  right-hand  side  of  Eq. 
(Al)  are  obtained  from  standard  tables"^  [aj  is  just  the  in¬ 
tegral  of  C„Hh)  over  the  lowest  300  meters  of  the  atmo¬ 
sphere  normalized  to  i j].  Table  1  shows  the  values 

of  hi  determined  by  Gauss-Legendre  quadrature  for  a 
four-screen  model  and  a  five-screen  model.  It  should  be 
noticed  that  the  lowest  phase  screen  is  always  placed  at 
an  altitude  of  150  m  to  ensure  sampling  of  boundary  layer 
turbulence.  (While  this  may  not  be  particularly  crucial  for 
an  investigation  of  stellar  scintillation,  it  nevertheless 
seemed  pinident  to  account  for  it  in  the  simulation.)  The 
strength  of  turbulence  for  the  various  phase  screens  may 
then  be  parameterized  by  a  set  of  Fried  parameters  J: 

/-va  ^  0A2Sk'^aiC,,Hhi).  (A2) 

Once  the  is  determined,  a  set  of  random  phase 

screens  is  generated  with  Fried’s  methodology.^^  Briefly, 
on  a  square  M  X  M  pixel  screen  we  generate  mean- 
zero,  delta-correlated,  Gaussian  random  variables, 
wijn^yrriy)  =  w{m},  with  a  variance  of  a^.  The  tw^o- 
dimensional  discrete  Fourier  transform  of  ii;(m),  W(  >/),  is 
then  multiplied  by  1 (except  at  j  77I  =  0,  where  W(  y) 
is  set  equal  to  zero)  to  produce  a  random  field  with  a  Kol¬ 
mogorov  spectrum.  Taking  the  inverse  Fourier  trans¬ 
form  of  yields  the  random  phase 

field  with  unconnected  tip  and  tilt  (to  be  discussed  below). 
The  variance  a“  is  chosen  so  that  the  structure  function 
associated  with  ^/^eChi)  matches  that  determined  by  the 
Fried  parameter  :  Uj  =  CM(Llr„  where  L  is 
the  length  of  the  phase  screen  and  C  =  0.1513128. 


Table  1.  Altitudes  of  Phase  Screens 


Screen  Number 

hi  for  4-Screens  fkm] 

hi  for  5-Screens  [km] 

1 

0.15 

0,15 

2 

2.52 

1.67 

3 

10.15 

6.80 

4 

17.78 

13.5 

5 

18.63 

In  generating  the  random  phase  field, 
ploy  a  screen  of  finite  size.  However,  strict  homogeneity 
of  the  random  field  requires  a  screen  of  infinite  dimen¬ 
sion,  implying  that  the  lowest-order  Fourier  components 
of  ^ 

similar  problem  exists  in  the  simulation  of  a  random  time 
series,  where  initiation  of  the  series  results  in  a  nonsta¬ 
tionary  random  process,^"^  Fried  corrects  this  error  by 
subtracting  the  incorrect  tip  and  tilt  from  and 

then  adding  a  statistically  appropriate  degree  of  tip  and 
tilt: 


^^(m,  hi)  ~  ^tte(^) 

-  (m 


M  \  j  M  \ 


(A3) 


Here,  6>(m,  hi)  is  the  random  phase  screen  employed  in 
the  computations,  I  is  a  unit  vector,  g  is  a  normal  random 
vector  (with  mean-zero  and  unit-variance  components),  y 
is  the  standard  deviation  of  the  statistically  appropriate 
tip  and  tilt  vector,  and  is  the  erroneous  tip  and  tilt  vec¬ 
tor: 


12 

M^iM^  - 


(M  +  1)  \ 

- ^ - 1) 

(A4) 


Fried  derives  the  statistically  appropriate  variance  of  tip 
and  tilt  from  the  phase  structure  function:  y^ 

=  2.610111/M(L/r,,,)^^^ 

Though  the  methodology  of  Fried  ensures  that  each 
phase  screen  can  approximate  an  appropriate  stochastic 
scalar  field,  Knepp^'*’  has  pointed  out  that  numerical  as¬ 
pects  of  the  simulation  problem  must  also  be  considered 
in  order  to  obtain  realistic  results.  From  the  form  of  the 
phase  structure  function  for  the  ith  screen,  we  have 
Z)a,^((5r)  =  where  Sr  is  the  phase 

screen’s  pixel  size:  dr  =  L/M.  Clearly,  if  the  screen  is 
to  accurately  represent  turbulent  refractive  index  fluctua¬ 
tions,  then  the  phase  variations  from  pixel  to  pixel  should 
be  less  than  90°.  We  achieve  this  constraint  by  requiring 
DoJ2Sr)  <  (77/2)^,  which  implies  an  upper  bound  on  5r 
that  is  determined  by  the  .  Since  is  always  larger 
than  where  we  get  a  criterion  for 

the  pixel  size:  5r<(r^/2)C 77^/27.54)^'^;  this  may  be  called 
the  structure-function  criterion  on  the  pixel  size.  Addi¬ 
tionally,  if  scintillation  in  the  saturation  regime  is  to  be 
accurately  simulated,  then  the  pixel  size  should  also  be 
less  than  the  Fresnel  length  associated  with  interscreen 
pi'opagation:  Sr  <  where  Azj.  is  the  smallest 

interscreen  distance;  this  may  be  called  the  interscreen- 
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Fresnel -length  criterion  on  the  pixel  size.  The  smaller  of 
these  two  criteria  places  an  upper  limit  on  the  pixel  size. 
In  the  present  simulations,  we  found  that  the  upper 
bound  on  dr  for  the  HA  profile  was  determined  by  the 
interscreen-Fresnel-length  criterion,  and  was  relatively 
insensitive  to  our  wind  speed  changes  in  turbulence. 
However,  for  the  HV-21  profile,  Sr  was  primarily  deter¬ 
mined  by  the  structure-function  criterion  and  did  change 
with  r  from  a  maximum  of  0.4  cm  to  a  minimum  of  0,3 
cm.  Note  that  as  a  general  rule  the  structure-function 
criterion  scales  like  whereas  the  interscreen-Fresnel- 
length  criterion  scales  like  (X/Ns)  ^^^. 

In  addition  to  an  upper  bound  on  the  pixel  size,  there  is 
also  a  lower  bound.  As  mentioned  above,  wave  propaga¬ 
tion  between  phase  screens  is  accomplished  by  Fourier 
transform  and  the  Fresnel  approximation.  Specifically, 
if  ij/ih)  is  the  two-dimensional  Fourier  transfoian  of  the 
field  at  /i,  then 


(A{/i2)  =  i/f(h))exp\ 


-  /12) 

2^ 


(A5) 


where  k  is  the  magnitude  of  spatial  frequency.  As  the 
values  of  k  are  discrete  in  the  simulation  (i.e.,  /c„i 
=  rmr/L,  in  =  1, 2,...M),  an  accurate  description  of 
propagation  with  use  of  the  Fourier  transformation  re¬ 
quires  that  the  difference  in  phase  associated  with  the 
two  largest  values  of  k  be  small  compared  with  ir. 


(M  - 

2^1^^ 


(A6) 


where  {h^  -  When  we  define  Lzi  as  the  larg¬ 

est  separation  between  phase  screens,  Eq.  (7)  results  in  a 
lower  bound  on  dy\  Sr  >  Note  that  to¬ 

gether,  the  upper  and  lower  bounds  on  Sr  constrain  the 
range  of  investigations  that  may  be  performed  for  a  given 
number  of  pixels,  M.  Under  conditions  where  the  upper 
bound  on  Sr  is  determined  by  the  iiiterscreen-Fresnel- 
length  ci'iterion,  M  must  be  greater  than  {2^Zi/^Zg), 
whereas  under  conditions  where  Sr^s  upper  bound  is  fixed 
by  the  structure-function  criterion,  M  must  be  greater 
than  (6.85XAz£,/r,,‘'^).  If  these  constraints  on  M  are  vio¬ 
lated,  then  there  is  no  value  of  Sr  that  satisfies  both  the 
upper  and  lower  bounds,  and  calculation  can  proceed  only 
with  an  increased  value  of  M. 

Finally,  because  of  the  finite  size  of  the  phase  screen, 
the  Fourier  transform  of  the  propagated  field  is  periodic, 
and  this  can  lead  to  spurious  interference  effects  near  the 
screens'  edges.  To  clarify  this  problem,  define 
,/ii)  as  a  field  element  at  the  plus  or  minus  trans¬ 
verse  side  of  an  upper-altitude  phase  screen.  Due  to  the 
turbulent  eddies  in  the  upper  screen,  the  plus  field  ele¬ 
ment  may  (for  example)  scatter  outside  the  footprint  of  a 
lower-altitude  screen,  whereas  the  minus  field  element 
may  stay  within  the  footprint.  However,  as  a  conse¬ 
quence  of  the  Fourier  transform’s  periodic  nature,  the 
plus  field  element  is  not  lost  in  the  calculation  but  reap¬ 
pears  at  the  opposite  edge  of  the  lower-altitude  screen  as 
a  spurious  field  element  [i.e., 

-  EspCp;,!,,.,  /i2),  while  hO  —  ECp^^ge*  ^2)]- 

The  total  field  at  the  minus  side  of  the  lower-altitude 
screen  is  the  sum  of  the  propagated  minus-side  field  ele¬ 


ment  plus  the  spurious  field  element:  £^(pedge»  ^2) 
-  ^2)-  Since  the  field  ele¬ 

ments  at  opposite  edges  of  a  screen  should  be  only  weakly 
correlated  (if  at  all),  this  effect  not  only  results  in  spuri¬ 
ous  partial  destructive  interference  at  the  screen  edges 
but  also  introduces  a  spurious  correlation  into  the  propa¬ 
gated  field  at  the  screen  edges.  To  guard  against  these 
edge  effects,  we  first  estimate  the  interscreen  scattering 
angle  for  propagation  from  screen  i  to  screen  i  +  1, 
6r-  where  is  defined  by 
=  1.  The  range  in  transverse  dimension  near  the 
i  +  1  screen’s  edge  that  is  corrupted  by  periodicity 
of  the  Fourier  transform  is  then  approximately 
(/ij  -  +  Adding  the  corrupted  ranges  for  all 

the  screens  provides  an  apodizing  radius,  r „  ,  at  the  tele¬ 
scope  aperture:  r^  ~  A/2  -  (hi  ~  hi  +  l)/ks^J  i .  As 

long  as  2ra  is  larger  than  the  telescope  aperture’s  diam¬ 
eter,  these  edge  effects  should  have  little  influence  on  the 
results. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs, 
specializing  in  advanced  military  space  systems.  The  Corporation's  Laboratory  Operations  supports  tiie 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staff’s  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure 
analysis,  solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared 
and  CCD  detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics, 
solid  state  laser  design,  micro-optics,  optical  communications,  and  fiber  optic  sensors; 
atomic  frequency  standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric 
propagation  and  beam  control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing 
and  evaluation,  battery  electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestmctive  evaluation,  component 
failure  analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures; 
launch  vehicle  fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics; 
chemical  and  electric  propulsion;  environmental  chemistry;  combustion  processes;  space 
environment  effects  on  materials,  hardening  and  vulnerability  assessment;  contamination, 
thermal  and  structural  control;  lubrication  and  surface  phenomena. 

Space  Science  Application  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
infrared  surveillance,  imaging,  remote  sensing,  and  hyperspectral  imaging;  effects  of  solar 
activity,  magnetic  storms  and  nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and 
magnetosphere;  effects  of  electromagnetic  and  particulate  radiations  on  space  systems; 
space  instrumentation,  design  fabrication  and  test;  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  tight  scattering,  state-specific  chemical 
reactions  and  radiative  signatures  of  missile  plumes. 

Center  for  Microtechnology:  Microelectromechanical  systems  (MEMS)  for  space 
applications;  assessment  of  microtechnology  space  applications;  laser  micromachining; 
laser-surface  physical  and  chemical  interactions;  micropropulsion;  micro-  and 
nanosatellite  mission  analysis;  intelligent  microinstruments  for  monitoring  space  and 
launch  system  environments. 

Office  of  Spectral  Applications:  Multi  spectral  and  hyperspectral  sensor  development; 
data  analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral 
imagery  to  defense,  civil  space,  commercial,  and  environmental  missions. 
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